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ABSTRACT 

Recent observations by the Large Area Telescope (LAT) onboard the Fermi satellite have revealed 
bright 7 -ray emission from middle-aged supernova remnants (SNRs) inside our Galaxy. These rem¬ 
nants, which also possess bright non-thermal radio shells, are often found to be interacting directly 
with surrounding gas clouds. We explore the non-thermal emission mechanism at these dynamically 
evolved SNRs by constructing a hydrodynamical model. Two scenarios of particle acceleration, either 
a re-acceleration of Galactic cosmic rays (GRs) or an efficient nonlinear diffusive shock acceleration 
(NLDSA) of particles injected from downstream, are considered. Using parameters inferred from ob¬ 
servations, our models are contrasted with the observed spectra of SNR W44. For the re-acceleration 
case, we predict a significant enhancement of radio and GeV emission as the SNR undergoes a tran¬ 
sition into the radiative phase. If sufficiently strong magnetic turbulence is present in the molecular 
cloud, the re-acceleration scenario can explain the observed broadband spectral properties. The 
NLDSA scenario also succeeds in explaining the 7 -ray spectrum but fails to reproduce the radio spec¬ 
tral index. Efficient NLDSA also results in a significant post-shock non-thermal pressure that limits 
the compression during cooling and prevents the formation of a prominent dense shell. Some other 
interesting differences between the two models in hydrodynamical behavior and resulting spectral 
features are illustrated. 

Subject headings: ISM: supernova remnants, shock waves, particle acceleration, 7 -ray 


1. INTRODUCTION 

Middle-aged SNRs have been receiving much attention 
recently thanks to the discovery of luminous 7 -ray emis- 
sion at their shells and sometimes close vicinities (e.g., 
Abdo et aL]|2009 20I0a|b|c Gastro fc Slane|[20l0 1 . The 
latest analysis of accumulated Fermi LAL data on the 


SNRs W44 and IG 443 by Ackermann et al. (2013) has 
clearly uncovered a kinematic cutoff at around 250 MeV 
in their 7 -ray spectra, a signature spectral characteris¬ 
tic that strongly supports the hadronic origin of the 7 - 
rays (i.e. production of neutral pions through proton- 
proton interactions and their subsequent decay into pho¬ 
ton pairs), and thus providing the long-awaited evidence 
for the acceleration of protons at SNR shocks. Earlier 
measurements b y AGILE have also shown evidences fo r 
such a signature (Tavani et al. ]Moi [Giuliani et al. J 2 ^. 
At the same tirne, however, the observed bright 7 -raiy 
emission raises new challenges to the conventional DSA 
theory of particle acceleration at astrophysical collision¬ 
less shocks. First of all, these dynamically evolved rem¬ 
nants have slow shocks and hence are not normally ex¬ 
pected to act as efficient particle accelerators like their 
younger cousins. This seems incongruent with the very 
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high 7 -ray luminosities inferred from observations. More¬ 
over, these 7 -rays that are concentrated in the GeV en¬ 
ergy band possess a characteristic spectral shape that 
deviates from a simple power-law with an exponential 
cutoff; rather, it points to the existence of a peculiar 
spectral break in the underlying proton distribution, typ¬ 
ically at a few to a few tens of GeV with a softening of 
spectral index by roughly one power from the lower to 
higher energy side of the break, the origin of which is 
still not clearly understood. Our current understanding 
and the major unknowns of the physics involved in parti¬ 
cle acceleration and non-thermal emission at a radiative 


shock is covered by a review by Bykov et al. 


(2013). 


Interestingly, the shocks of middle-aged biNRs are usu¬ 
ally found to be propagating through high-density en¬ 
vironments such as adjacent massive molecular clouds, 
evidenced by observational features like 1720 MHz OH 
masers and various optical and infrared emission lines. 
Signatures like strong forbidden line emission around the 
shocks suggest that the shocks have undergone transition 
into the radiative phaselj It is thus natural to ques¬ 
tion whether radiative shocks in molecular clouds can 
somehow manage to generate bright non-thermal radio 
and 7 -ray emission with the observed properties, despite 
their low velocities compared to non-radiative shocks at 
younger SNRs. Several theoretical models have already 
been presented in the literature even before the matu- 
rity o f 7 -ray astronomy. For example, [van der Laanj 
(19621 showed that compression of cosmic ray electrons 
and rhagnetic field in an isothermal shock could account 
for the radio synchrotron emission of older shell-type su- 


® A broader overview on the rich observational properties of 
these objects can be found, e.g., in the recent review article by 
[Slane et al.||^14||. 
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pernova remn ants such as the Cygnus Loop.|Blandford fc] 
'Cowie| (119821) extended this work, and |Cox et al.| p999|) 
applied th e idea to W44. With th e advent ot 7 -ray ob- 
servations, Uchiyama et al. (2010| (hereafter ‘UlO’) ad¬ 
vanced this scenario and constructed a phenomenologi¬ 
cal model using an analytical approach. They concluded 
that fast radiative J-shocks are able to pick up and re¬ 
accelerate pre-existing CRs. If these re-accelerated CRs 
are subsequently boosted further to higher energy density 
by compression inside a rapidly cooling and contracting 
shell behind the radiative shock, they can simultaneously 
produce sufficiently bright GeV 7 -ray and radio syn- 
chrotron emission to expla in observations. Most recently, 
Tang & Chevalier ( 2014b|) attempted to construct a time- 
dependent description of DSA at slow shocks based on 
a simple an alytic hydrodynamic mod el for shock-cloud 
interactions (|Tang & Chevalier||2014a| and a parametric 
model for the CR diffusion coefficient, and suggested that 
a time-dependent test-particle solution of re-acceleration 
under a Kraichnan-like diffusion coefficient {D{p) oc 
can explain the observed 7 -ray spectra of SNR IC 443 
and W44. 

In this paper, we investigate the immediate locality 
of one of these cloud shocks and explore its broadband 
non-thermal emission mechanisms using a fully time- 
dependent hydrodynamical simulation, self-consistently 
coupled to an explicit treatment of DSA (either re¬ 
acceleration of pre-existing CRs or direct acceleration 
of the downstream thermal particles) and other impor¬ 
tant physical processes including those associated with 
a radiative shock. We conclude that, despite minor dif¬ 
ferences, our results broadly confirm the assertions by 
UlO that a re-acceleration model can well explain the 
observed broadband emission properties. A direct ac¬ 
celeration model on the other hand suffers from several 
difficulties that we will explain in detail. In Section 
we introduce the essential physical components of our 
model and explain their relevance to the problem. Ma¬ 
jor assumptions made are also stated and elaborated. 
We then present our results on the hydrodynamical evo¬ 
lution and non-thermal emission calculation in Section 
followed by discussions and conclusion in Section]^ 

2. MODEL 

In this section, we review the physics and numerical 
setup of our model. The numerical calculations are per¬ 
formed usi ng the CR-hydro-NEI simulation code (Lee, 
et al.|2012 ). In our model, we consider a picture in which 
the progeiiitor star undergoes core-collapse and explodes 
into a ten uous wind cavity sur rounded by a molecular 
cloud (e.g. Koo & Heiles 1995). The blastwave created 
by the expanding ejecta propagates into the cavity with 
high speed and eventually hits the interface with the 
surrounding dense medium (probably part of the parent 
molecular cloud of the progenitor star), typically at a ra - 
dius ~ 10 pc from the explosion center (|Chevalier|1999 ). 
It then penetrates into the cloud (we will call it a ‘cloud 
shock’ hereafter) and keeps pushing through the dense 
medium, and eventually becomes radiative. The molec¬ 
ular cloud is approximated as a uniform dense medium 
with a number density no 100 cm“^. While the real 
situation is expected to be more complex with three- 
dimensional details such as dumpiness (i.e., density in- 
homogeneity) of the molecular cloud (e.g., Rho et al. 


1994), we believe our spherically symmetric model can 


capture the essential physical phenomena and provide a 
fundamental understanding of such systems. This step 
is an essential precursor to more complicated and time- 
consuming three-dimensional hydrodynamical or MHD 
simulations that self-consistently include DSA and other 
important physical processes. 

We will trace the time evolution of the hydrodynam¬ 
ics of the cloud shock starting from the time of pene¬ 
tration, passing the point when the shock has deceler¬ 
ated to about 200 km s“^ and become fully radiative. 
Physics including particle acceleration through the dif¬ 
fusive shock acceleration (DSA) mechanism and micro¬ 
physical processes such as ionization and recombination, 
thermal conduction, radiative cooling and photoioniza¬ 
tion heating are calculated simultaneously. At the same 
time, charge exchange can be an important effect in as- 
trophysical shocks. Concerning the physical conditions 
we are investigating here, the balance between ionization 
and recombination can be very sensitive to temperature 
and density, with additional complications due to heat¬ 
ing terms such as photoionization. Unfortunately, our 
model does not currently include an accurate treatment 
for the effects of charge exchange. We recognize that 
charge exchange could have an important effect on our 
results, and we will address those effects in a subsequent 
paper. This work also does not concern the reverse shock 
and the ejecta emission which are relevant to the thermal 
X-ray emission. Such discussions are reserved to ano ther 
line of studie s we are carrying out in para llel, such as |Lee| 
et al.l (|2014|) and iPatnaude et al.l (120151) , which will be 


(j2014|) and Patnaude et al. (|2015|), which will 
a with the current work in the future for a more 


mergeo 

complete description of middle-aged SNRs. 

The typical timescale at which the shock becomes ra¬ 
diative for a supernova (SN) explosion of energy U 51 = 
A'Sn/( 10 ®^ erg) inside a uniform surroundin g medium of 
density no is given by iBlondin et al. (1998): 


ttr ~ 2.9 X lO'^Ug^^^UQ yr 


( 1 ) 


By replacing Usn with the total energy of the driving 
blastwave in the cavity and no with the cloud density, 
this quantity is found to be a convenient time unit for 
measuring the radiative transition of the transmitted 
cloud shock as well, provided that effects of efficient par¬ 
ticle acceleration and other sources of energy loss which 
can further speed up the deceleration of the cloud shock 
are not significant. Hence, we use it as the basic time 
unit for our models throughout this work. Rapid cool¬ 
ing occurs at about t > Ur at which a cold, condensed 
dense shell should form behind the radiative shock where 
it is pushed outward by the hotter interior. Compression 
of the gas and the transverse component of the R-field 
lines is expected inside the contracting cool shell. The 
same should also occur with the trapped (re-)accelerated 
CR particles inside the shell which not only increase in 
density but are also energized by the compression, i.e.. 


Ap = (s^/^ - l)p , (2) 


where Ap is the momentum increase of a CR particle 
experiencing the compression and s 1 is the compres¬ 
sion ratio in the cool shell. We assume that the R-field 
is frozen into the plasma behind the shock. 

The enthalpy density id in a Lagrangian gas element 
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can be written as: 

TT ^9 r> I TCR I -^_L I ^ 2 (o\ 

H = —^Pth + - tPcr + ^ + 7;Pv , (3) 

7g - 1 7 cr - 1 47r 2 

where p is the mass density, v is the flow speed, i3j_ is the 
B-field component perpendicular to the shock normal, 
Pth and PcR are the thermal and CR pressures, jg and 
7CR are the ratio of specific heats for ideal gas and CRs 
respectively. The evolution of enthalpy per unit mass 
h = H/p in the gas e lement is f ollowed using the energy 
equation below (e.g., Cox||1972), 


dh 

dt 


—UeriH^ix) + nr(a;) + kV^T 


( 4 ) 


where nf.^ uh and n are the electron, hydrogen and total 
number densities respectively, and A, T and k are the 
cooling coefficient, heating coefficient and thermal con¬ 
ductivity respectively. The treatment of each of these 
components will be described below in more detail. 

In the immediate post-shock region, it has been sug¬ 
gested by observations that the electron-to-proton tem¬ 
perature ratio is small for shocks faster than 1000 km s“^, 
but close to 1 for slowe r shocks such as those we are 
interested in here ( e.g., iGhavamian et al.||20011 [2007 


Rakowski et al. 2008 Helder et al.|2010 ). We hence adopt 


an instantaneous post-shock equilibration of tempera¬ 
tures among ions and electrons. The elemental abun¬ 
dances in the dense cloud are taken from the observations 
and models of bright at omic lines from SNR molecular 
shocks in the IR band by Reach & Rho (20001. All molec¬ 
ular species are expected to be fully dissociated by the 
fast cloud shock. 

2 .1. Non-equilibrium Ionization, Cooling and Heating 

The fully time-dependent ionization and recombina¬ 
tion of 12 chemical elements including H, He, C, N, O, 
Mg, Ne, Si, S, Ar, Ca and Fe are followed behind the 
shock. The number densities of a total of 152 ion states 
are calculated in each gas element and time step in order 
to estimate effects like radiative cooling and photoioniza¬ 
tion heating self-consistently. The details of our treat- 
m ent of non-equilibrium ioniza tion (NEI) can be found 


Patnaude et al. (2009 2010). 


3y tracing the ionization states and number densi¬ 
ties of all ion species in the post-shock plasma through 
the NEI calculation, we can compute their contributions 
to radiative cooling through various fine and continuum 
emissions at each time step. This approach is more accu¬ 
rate than applying a power-law cooling curve which only 
depends on the temperature. A selection of 12 relatively 
strong optical and infrared collisional excited lines (in¬ 
cluding [C II] 156/im, [N II] AA6548, 6584, [O I] AA6300, 
6364, [O II] AA3727, 3729, [O III] AA4959, 5007, [O IV] 
25.8/rm, [Ne II] [Si II] 35Atm, [S II] AA6717, 6727 

and [Ee II] 26/rm) are also computed. In the model con¬ 
sidered in this work, the cooling timescale is typically 
much smaller than the dynamical timescale (sound cross¬ 
ing time) in the post-shock region of rapid recombination 
and radiative cooling, hence the cooling there can be ap¬ 
proximated essentially as an isochoric process. This re¬ 
gion is then subsequently pushed upon and compressed 
by the faster and hotter interior to form a dense shell un¬ 
til a mechanical equilibrium is approached. This equilib¬ 


rium in the dense cool shell is predominantly supported 
by the non-thermal (i.e., magnetic and cosmic-ray) pres¬ 
sures, as will be discussed further later in this paper. 
During the transition to the pressure-driven phase, large 
amplitude oscillations of the radiative cloud shock due to 
thermal instability are expected in an over-stable fashion 
until the shock velocity becomes lower than ^ 120 km s“^ 


1997 


Hlondin et al.||1998]), as we will also clearly show 
rde 


Chevalier & Imamura 1982 Kimoto & Chernoff 


below) We do not consider cooling and heating effects 
from molecular chemistry like molecule reformation and 
IR line emission below a few 100 K. 

Absorption of and photoionization by locally generated 
UV photons can heat up the gas by the production of 
photoelectrons. For simplicity, we assume a temperature 
of To ~ 10^ K in the precursor due to photoionization 
heating after the shock becomes radiative. We also do 
not treat the time-dependent radiative transfer in the 
post-shock gas. We calculate the heating rate r(x) by 
estimating the fluxes and optical depths of photoionizing 
photons including the continuum emission and important 
emission lines, such as He I, He H, O HI-V etc, in each 
gas parcel. The absorpti on cross sections are calculated 
using the fitting recipe of Verner et al. (|1996|). Follow ing 
a method similar to that in Cnat & Sternberg (2009), we 
calculate in a steady-state limit the absorption of these 
self-generated photons by each chemical element in the 
gas cells further downstream and the resulting produc¬ 
tion of photoelectrons whose residual energies heat up 
the plasma. 

Thermal conduction can be important at places with 
large spatial temperature gradients, such as the rapidly 
cooling gas behind a radiative shock propagating in a 
molecular cloud. We treat thermal conduction by pa¬ 
rameterizing the conductivity relative to the Spitzer con¬ 
ductivity KSpitzer as 


/ _ 7 ^ 5/2 \ 

’ — /cond ^Spitzer — /cond [1.84 X 10 X J 


( 5 ) 


where InAc « 37 is the Coulomb logarithm. In our 
model, we set th e parameter fmnri = 0.3 for a collisionless 


plasm a (see e.g. Narayan & Medvedev 2001 Bale et al. 
2013), but our results are not sensitive to a moderate 
change to this parameter. 


2 . 2 . 


Non-thermal Emission from a Cloud Shock 

Currently, there still remain controversies on the true 
origin of the center-filled thermal X-ray emission from 
middle-aged, mixed-morphology SNRs. For example. 


ogy (e.g., [White & Long 

1991 

whether they 

are dominated t 

matter (e.g., 

Zhou et al. 

2014 

Uchida et al.jj 

2 U 12p which most 


Cui & Cox 1992), and 


probably have to be dis- 


going debate on the physical explan ation for the appar¬ 
ent “over-ionized” plasma state (e.g. , Itoh fc Masai|1989 


Kawasaki et al.|2002[ |Zhou fc Ciie^|2011| Shimizu et al.l 

2012 |), evidenced by the detection of radiative recombina- 

tion continuum components and accompanying anoma- 
lous ratios of H-like and He-like ionization states (e.g. 


Kawasaki et al.|2002 Ozawa et al. 12009) [Yamaguchi et ah 

2009| ). The situatioii tor the non-therrnal emission is ar- 
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guably more clearcut. The major contributions to the 
bright radio continuum and 7 -ray emissions, both pre - 
sumably coming from the shell (e.g. Abdo et al.|| 2010 b ), 
are synchrotron radiation emitted by primary (i.e., di- 
rect acceleration of thermal e“ or re-acceleration of pre¬ 
existing CR e“ and e"*") and secondary (i.e., e~ and e+ 
from the decay of changed pions produced by hadronic 
interactions) leptons that gyrate around the post-shock 
B-field, and energetic 7 -ray photons from the decay of 
neutral pions produced by the interaction of accelerated 
protons/ions with the dense gas around the shock, re¬ 
spectively. We will adopt this general picture in our 
models for the non-thermal emission. 


2.3. Particle acceleration and re-acceleration 

There are two possible scenarios that we can consider: 
( 1 ) re-acceleration of pre-existing high-energy particles 
by the cloud shock; ( 2 ) efficient acceleration of particles 
through NLDSA via injection from the thermal pool be¬ 
hind the subshock. 

The re-acceleration scenario involves the shock accel¬ 
eration of pre-existing high-energy particles wandering 
around the source, such as Galactic CRs, and escaped 
CRs from nearby younger SNRs. In this case, most of the 
pre-accelerated “seed” particles already have adequate 
momenta to avoid the difficulties of thermal injection due 
to Coulomb loss and long acceleration timescales. One 
difficulty can be the relatively low densities of the seed 
particles compared to the observed high luminosities of 
the radio and GeV emission. We will however show that 
this hurdle can be overcome if the shock undergoes a 
transition into the radiative phase. 

In this model, following UlO, we assume the seed 
particles to be the Galactic CR protons and elec- 
trons-bpositrons (hereafter we refer to them simply as 
electrons) and adopt spectra in momentum space in the 
following forms, 

^p,seed(?3) = 47rJp/3 Pq , Po b> 0.31 , 

ne.seed(p) = 47r Jg (1 -b Pg) > PO > 0-02 , 


where po = p/(GeVc“^), /3 = Vi/c, Vi is the particle 
velocity. Here nseed(p)dp is the particle number den¬ 
sities in the interval p p dp which can be rewrit¬ 
ten as ‘in’p^fseed{p)dp in terms of the phase space dis¬ 
tribution function /seed(p)- The normalization factors 


adopted are Jp = 1.9 cm 


and Jp = 


0.02 cm ^s ^sr ^GeV so the total energy densities 

t-3 


Jp 

‘■s^'sr” 

of the seed protons and electrons are Sp « 0.81 eV cm” 
and Ee « 4.8 x 10“^ eV cm“^ respectively. With such, 
we can calculate the spectra of the re-accelerated parti- 
cl es using an i terat ive semi-analytic approach described 


Blasi (20041 and Lee et al. (2012), as follows 


fi{p) = 


iSu 


5totC/(p) - 


X exp 


X exp 


r dp"{ 

T/ 777) A.seedlP ) 

-'Pmin ^ y 

fP dp' 3StotU{p') 

Ip" p' 5'totG(p')-l 


P 


Pmax,2 


(7) 


where i represents e”/"*" or p, U{p) and S'tot are the di¬ 
mensionless gas flow velocity in the shock rest fram^ 
and the compression ratio experienced by the streaming 
GRs in the shock precursor respectively, Pmin is the cor¬ 
responding low-energy cutoff of the seed spectra specified 
in Equation (|^, and Pmax is the maximum momentum. 
The index a specifies the rollover shape near Pmax and 
we adopt a = 1.5 for o ur models. For de tailed definitions 
of these quantities, see Lee et al. (2012 1 . We do not con¬ 
sider the acceleration of heavy ions or dust grains here. 

For the NLDSA scenario, we expect that DSA at 
the cloud shock with a velocity Ugk ~ 100 km s“^ 
and density ng > 100 cm“^ may face difficulty if the 
so-called ‘thermal leakage’ mechanism is considered for 
the injection of suprathermal particles into the DSA 
process. A first possible difficulty is the fast energy 
loss of the injected particles through Goulomb inter¬ 
actions and ionization which competes against their 
early acceleration and hence hinders their path towards 
higher energies. This constraint is irrelevant if the ratio 
tacc/min(tcoui,iion) < 1 , where face, ^Coui and tio„ are 
the acceleration time-scale, Goulomb loss time-scale and 
ionization loss time-scale respectively at around the in¬ 
jection energy. For a temperature ^ 10'^ K and H-field 
strength ~ 10 pG, this condition can be expressed ap- 
proximately as (us k /100 km s“^)^/(no/l cm“^) 10 “^ 

( Drury et al.|l99(3 ), which does hold for ng ~ 100 cm s ^ 
and Vsk ~ 100 km s“^. Furthermore, in situations where 
the cloud is highly clumpy, the shock can be propagating 
in an inter-clump medium (e.g., Tang & Chevalier 2014aI 
which has a much lower average density than the cloud 
cores and hence poses less problem from Coulomb and 
ionization losses. 

However, even if Coulomb and ionization losses are not 
problematic, the slow shock speed still implies a long ac¬ 
celeration time required for boosting the particles from 
the injection energy up to 7 -ray emitting energies, typ¬ 
ically > 100 GeV. One possibility to overcome this is 
to consider a situation with significant amplification of 
the magn etic turbulence, su ch as by CR-driven instabili¬ 
ties (e.g., Bykov et al.||2014), to shorten the acceleration 
time. In this paper, we are going to discuss such a case 
with an efficient NLDSA of particles injected from the 
thermal pool, which can produce strongly amplified B- 
field through CR steaming instability in the shock pre¬ 
cursor. Our formulation for NLDS A and its feedbac ks to 
the hydrodynamics can be found in Lee et al. (2012 1 and 
references therein. 


2.3.1. Diffusion Coefficient and Momentum Break 

It is straightforward to confirm that a Kolmogorov or 
Kraichnan magnetic turbulence spectrum typically asso¬ 
ciated with the interstellar medium (ISM) or the insides 
of giant molecular clouds {D ~ 10^^ to 10^® cm^ s“^ at 
10 GeVc“^) is unable to support acceleration of particles 

^ P{p) 1) normalized to the shock speed, and is a func¬ 
tion of X but with a change of variable to p using the momentum- 
dependent diffusion length of a CR particle in the shock precursor. 
It represents the smoothed shock structure when CR back-pressure 
strongly modifies the incoming flow speed profile. For example, 
U(p) Si 1 for all p if shock modification from CR pressure is unim¬ 
portant, and [/(p) < 1 when DSA is efficient due to a deceleration 
of the incoming flow by the back-streaming CRs in the shock rest 
frame. 
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(a) 

Fig. 1.— (a) Spectra of re-accelerated Galactic CR at a steady cloud shock of velocities Usk < 120 km s“^. Here, the ambient magnetic 
field B n = 25 fiG and Tq = 10^ K in the shock precursor. The precursor ion fraction /;□„ is estimated according to |Hollenbach &: McKee| 
(|1989|l which changes the break momentum of fho re-accelerated particles. The dashed lines show the ambient spectra ot the GKs before 
re-acceleration. The upper thick and lower thin lines correspond to protons and electrons respectively, (b) Variation of esacpe-limited 
maximum momentum pmax, break momentum precursor ion fraction /ion with shock velocity Ugk in correspondence with the 

spectra shown in panel (a). 


to energies that can generate substantial sub-TeV pho¬ 
tons at a slow cloud shock in a dense medium. It has 
been shown that stronger turbulence which reduces the 
diffusion coefficient by a factor of ^ 100 relative to the 
Galactic mean is necessary to reprod uce the observed 7 - 


ray fluxes of middle-aged SNRs (e.g., Fujita et al. 2009 


Ohira et al. 2011 Yan et al. 2012). Observationally, there 


is evidence for such turbulent shock precursors from re¬ 
cent Ho obs ervations of SNR shocks interacting with gas 
clouds (e.g., Lee et al!]|2010 Medina et al.||2014 |. 

One possibility tor the presence ot strong magnetic tur¬ 
bulence in the molecular clouds is an efficient, non-linear 
CR acceleration through DSA at the cloud shock which 
can amplify the turbulent component of the upstr eam B- 


tum pbr above which wave-particle scatterings become 
weaker and less efficient. 

At around pbr, the CR diffusion in the precursor 
changes from Bohm to a regime in which the CRs scatter 
only with small-scale waves and can be described by a 
diffusion coefficient D{p) oc In this light, once the 
shock velocity has dropped to Usk < 120 km s“^, we 
adopt a diffusion coefficient for the CRs with a smooth 
break at pbr with the following form: 


D{x,p) = rjsDBix.p) ( l-f — 

Pbr 


( 8 ) 


field throug h various CR-driven instabilities (e.g., Bykov 
et al.|[2M4|, which applies to our NLDSA scenario. An¬ 


other possibility is that efficient acceleration occurred in 
the past when the SNR shock was propagating at high 
velocities in the wind cavity, from which high-energy 
CRs that escaped upstream from the shock can pene¬ 
trate through the surrounding cloud and generate strong 
magne tic turbulence th ere through the streaming insta¬ 
bility (Yan et al.||2012 1. 

In either scenario, we assume that the B-field in the 
molecular cloud is highly turbulent and fully tangled to 
support Bohm-like diffusion of CRs up to a certain par¬ 
ticle momentum above which the ion-neutral damping 
effect in a partially neutral medium takes place. At the 
cloud shock, the upstream medium can be expected to 
be partially ionized. Even after the shock has become 
radiative and the UV light produced in the downstream 
can create a photoionization precursor, the ionization is 
partial and neutral particles exist for Usk ^ 120 km s“^. 
In such a case, ion-neutral damping of the magnetic tur¬ 
bulence in the shock precursor will occur and hamper CR 
acceleration at the highest momenta. More specifically, 
the pitch-angle scatterings of the CRs with the turbulent 
B-field will experience a transition at a certain momen- 


where Dg is the Bohm diffusion coefficient and pg ~ 
1. For a radiative shock with higher velocities, full 
ionization is assumed in the precursor and D(x,p) = 
VbDb{x,p) is used. 

We can est i mate the break momentum by following 
Maikov et al. (2011) and references therein, i.e.. 


Ph. 


= 10 ( 


Tn 


y 

Vio^k; \pG 


-1 - 1/2 
n„ n, rupc 


(9) 


where and rii are the neutral and ion densities. The 
ion fraction /ion = nijirii -I- n„) in the precursor of the 
cloud shock can be estimated approximately from Vsk us¬ 
ing the results of Hollenbach & McKee (19891. A spectral 
break for the CR spectra can be expected if pbr < Pmax, 
and we can impose the break b y multiplying a factor of 
(ip/i Ohr)~^ to fi(v) for jo > pbr- Draine & McKee (1993) 


and Drury et al. (1996) also estimated the maximum en¬ 
ergy attaInaBIeD 3 rtEieaccelerating CRs when ion-neutral 
damping limits DSA. The maximum energies they de¬ 
rived have the same ion fraction dependence as phr in 
Equation 

The ion fraction /ion is estimated approx imately us¬ 
ing the st e ady-fl ow equilibrium value from |Hollenbach| 
& McKee (1989) in our model. We note however that 
non-steady-state effects can bring about important mod- 
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ifications. For example, | Cox (1972) has pointed out that 
an enormous burst of EU V photons when the cool shell 
forms can photoionize the precursor, and the recombina¬ 
tion timescale is longer than the time for further slowing 
of the shell. This can lead to a more highly ionized gas 
compared to the equilibrium values at < 120 km s“^. 
Also, /ion can change if non-linear feedback of PSA t o 
the shock structure is considered (iBykov e t aI1| 2013|), 
and the formulae in Hollenbach & iVlcKee| ( 1989P must 
be modified; shock rnoditication by the Cll pressure re¬ 
sults in an increase of the total compression ratio and 
decrease of the subshock compression ratio, which leads 
to a reduction of the post-shock temperature relative to 
the standard value for unmodified shocks. This hastens 
the transition of the shock to the radiative phase, and 
the ion fraction can start to increase due to photoioniza¬ 
tion at a higher shock velocity. In our models, however, 
DSA does not enter the nonlinear regime for the sce¬ 
nario of Galactic CR re-acceleration. And for the case of 
NLDSA of thermally injected particles, the evolution of 
the ion fraction in the precursor is expected to be highly 
non-trivial, especially when the time-dependent NLDSA 
feedback to the shock hydrodynamics is accounted for. 
In this work, we neglect this nonlinear effect on the ion 
fraction for simplicity. 

2.3.2. Spectrum of Accelerated Particles 

In the case of re-acceleration, the maximum momen¬ 
tum Pmax of the CR protons at the cloud shock is 
most likely determined by their escape far upstream. 
We model their escape through setting a free-escape- 
boundary (FEB) at 10% of the current shock radius 
ahead of the shock, so that CRs with diffusion lengths 
longer than the FEB in the precursor are considered to 
have escaped. 

For particle acceleration via injection from the ther¬ 
mal pool, Pmax is additionally limited by the finite ac¬ 
celeration time due to the slow shock speed in a molecu¬ 
lar cloud; in fact, significant magnetic field amplification 
through CR-driven instabilities in the precursor is neces¬ 
sary to shorten the acceleration time so that the injected 
thermal particles can reach sub-TeV energies in a rea¬ 
sonable time, therefore efficient NLDSA is required. 

When ion-neutral damping in the shock precursor is 
important, either before the shock becomes radiative 
or as the radiative shock has decelerated to a velocity 
<120 km s“^, Pmax is limited due to the transition of 
D{x,p) to a fast regime above the break momentum (i.e., 
waves with frequencies w = kvA smaller than the ion- 
neutral collision frequency izi„ become strongly damped, 
with k and va being the wave number and Alfven ve¬ 
locity). The Pmax of the CR electrons can be further 
limited by energy loss processes such as synchrotron and 
inverse-Compton losses if important. For more details on 
the calculation of Pmax, see the descriptions in|Lee et al.| 

([M^. 

As an example. Figure]^ (left) shows the spectra of 
re-accelerated Galactic CR protons and electrons at a 
steady-state cloud shock with different Wsk in a dense 
medium with ambient proton density uq = 200 cm“^, 
magnetic field Bq = 25 /rG and temperature Tq = 10^ K. 
The pre-shock ion frac tion /inn of the incom i ng ga s is 
calculated according to Hollenbach & McKee ( 1989|) as 
described above. The cloud shock radius trorn the ex¬ 


plosion center is set at 10 pc here, but the result is not 
sensitive to a fractional change of this value. 

Figure (right) shows the corresponding variation of 
Pmax, Pbr and /ion as a function of Usk- Here Pmax is 
limited by escape and is determined by equating the CR 
diffusion length to the distance between the FEB and 
the shock front through Equation]^ The maximum mo¬ 
mentum of electrons (not shown) is very close to that 
of the protons since magnetic field amplification due to 
self-generation of magnetic turbulence by the CR stream¬ 
ing instability in the shock precursor is inefficient in the 
case of pure re-acceleration of Galactic CRs, so that en¬ 
ergy loss of the re-accelerating electrons via synchrotron 
radiation is slow in the immediate downstream region. 
However, synchrotron losses can be much faster further 
downstream when a highly compressed H-field is present 
in the radiatively cooled dense shell formed behind the 
cloud shock. 

In this model, the precursor of a radiative shock with 
Vsk > 120 km s“^ approaches full ionization by photoion¬ 
ization, and wave damping through ion-neutral collisions 
is ineffective, so no break in the diffusion coefficient and 
the CR spectra appears. At lowest velocities, /ion "C 1 
and the ion-neutral collision frequency in the precur¬ 
sor is small, hence wave damping becomes less efficient 
and the break momentum pbr can rise again according 
to Equation Therefore, Pmax does not scale trivially 
as a simple power of Vsk in this velocity regime. For pa¬ 
rameters adopted in this example, a minimum for pbr 
is reached at v^k ~ 100 km s“^ where there is a com¬ 
parable concentration of ions and neutrals in the pre¬ 
shock gas. However, we note that this semi-analytic re¬ 
sult is obtained by considering steady shocks only, and 
can be modified if detailed wave-particle and particle- 
particle interactions are self-consistently followed for a 
time-evolving radiative shock, which is beyond the scope 
of this paper. 


2.4. Shock Hydrodynamics in a Dense Medium 

Using the re-acceleration model as an example, we can 
first look at the general hydrodynamical behavior of a 
cloud shock just before and after its transition into the 
radiative phase. Here and throughout this work, we 
adopt a distance dsNR = 2.9 kpc to the SNR, a final 
SNR radius i?sNR = 12.5 pc, an ambient gas density 
no = 200 cm“^, magnetic field Bq = 25 ^G and pre¬ 
shock temperature Tg = 10"^ K. The pre-shock density 
chosen is typical of the molecular shocks at SNR W44 
which are most probably responsibl e for producing the 


bright radio synchr otron filaments (|Reach et al. 2005 
Yoshiike et al. 20131. 


'Lhe evolution of the shock hydrodynamics in the dense 
uniform medium is followed up to a time Time snap¬ 
shots of the profiles of hydrodynamic variables includ¬ 
ing total gas density n, temperature T, magnetic field 
strength B and gas velocity v are shown in Figure until 
the shock has just started to become radiative. When the 
shock decelerates quickly with time in the dense medium 
and the post-shock gas temperature drops to around a 
few 10® K, radiative cooling becomes important. An en- 


® Note that this ‘cloud shock time’ tc is counted from the pen- 
etration of the blastwave into the cloud, so it is smaller than the 
actual age of the SNR. 
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(a) Time Snapshots of Hydrodynamic Profiles 



(b) Evolution of Cloud Shock Velocity 

Fig. 2. — (a) Time snapshots of profiles of hydrodynamic vari¬ 
ables up to the stage when the shock has just become radiative for 
the re-acceleration model. Prom top to bottom: gas temperature, 
gas (proton) density, H-field strength and gas velocity. Radius is 
in units of the FS radius. The colorbar depicts the evolution age 
of the cloud shock in units of characteristic transition time ttr. (b) 
Cloud shock velocity Ugi; as a function of time for the re-acceleration 
scenario. 


hancement of density and i3-field with time due to rapid 
compression can be observed. A cool dense shell forms 
behind the shock which is pushed forward by a faster and 
hotter interior heated up earlier by a stronger shock in 
the past. In this phase, the dynamics approaches that of 
a pressure-driven blastwave. The cool shell is eventually 
dominated by and supported against further collapse by 
non-thermal pressures, as is known to be the case at ra- 
diative filaments in som e remnants such as Cygnus Loop 
(Raymond et al.||1988). The formation of the dense cool 
shell can happen at a time somewhat later than 1 ttr 
due to the inclusion of thermal conduction which redis¬ 
tributes heat to slow down the cooling process. After 
the beginning of the transition, the cloud shock shows 
large oscillations in velocity for a few cycles due to ther¬ 
mal ins tabil ity (Figure^ bottom panel), as mentioned in 
Section [2T| The instarality subsides and the oscillations 
are damped after the mean shock speed becomes lower 
than ~ 120 km s“^. 


3. RESULTS 

In this Section, we consider a system where the SN 
explosion occurred in a wind cavity in the past and the 
blast wave eventually hit a surrounding molecular cloud. 
We follow the dynamics of the transmitted cloud shock 
and calculate the associated radio and 7 -ray emission 
under two scenarios of (i) Galactic CR re-acceleration, 
and (ii) NLDSA of thermally injected particles. To see if 
they can reproduce the generally observed characteristics 
of GeV-bright middle-aged SNRs, a comparison of the 
models with broadband spectral data from SNR W44 is 
presented. 

3.1. Case of Re-acceleration 
3.1.1. Broadband Non-thermal SED 

Figure shows the spectral evolution of the radio and 
7 -ray emission predicted by our model under the Galac¬ 
tic CR re-acceleration scenario. The corresponding evo¬ 
lution of integrated energy fluxes is also shown. For clar¬ 
ity, the contribution from non-thermal bremsstrahlung 
to 7 -rays is not included in the plot since it is relatively 
unimportant compared to the 7 r*^-decay component. The 
time tc is counted in units of ttr for generality. The av¬ 
erage Usk right after the penetration is <400 km s“^. 
At tc < 1.2 ttr, both the radio continuum and 7 -ray 
flux increase gradually with time, as more re-accelerated 
CRs and secondary from pion decay accumulate 

downstream. The spectral shapes do not change signifi¬ 
cantly due to the gradual evolution of the shock. Then, 
a prominent transition occurs at tc « 1.2 ttr when the 
shock has decelerated to a speed < 200 km s“^ whose 
dynamics start to be governed by radiative cooling, just 
as we have seen in Figure]^ The compression of the B- 
field and gas density as well as the re-energization of the 
re-accelerated GRs and secondary particles in the rapidly 
cooling and contracting gas shell boost the non-thermal 
fluxes both in radio and 7 -ray by a factor of a few 100 . 
The spectral shape of the 7 -rays also changes as the shock 
becomes slower than 120 km s“^ and experiences a par¬ 
tially ionized precursor so that the particle acceleration is 
hampered by the ion-neutral damping effect we discussed 
above. 

Eventually, as the gas cools down and compresses, the 
non-thermal pressure becomes higher than the thermal 
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(a) Evolution of Broadband Spectrum 



(b) Evolution of Integrated Flux 

Fig. 3.— (a) Time evolution of the radio synchrotron (left panel) 
and TT^-decay 7 -ray (right panel) spectra up to tc = 3 ttr of our 
model under the Galactic CR re-acceleration scenario. Dashed 
lines in the left panel show the contribution from secondary e~ 
and e"*" to the synchrotron emission, while the solid lines show 
the total. Again, the colorbar depicts the evolution age of the 
cloud shock in units of ttr- Observed data of SNR W44 from radio 
telescopes (compilation by ICastelletti et al.| 2007l ISun et al.|2011| 
black points ), Planck Galactic blMH survey ([Pianck tJoilaboration| 
|et al. | 2014| grey points) and 7 -ray instrunients includi ng t'ermi 
LiAi' ^^gri ^nn et ah |2013| black points) and AGILE (|Cardillo| 
|et al.pUi4| green points) are overlaid. The radio data points are 
scaled by a factor of 0.5 following UlO. An overall normalization 
factor of 0.2 is applied to the model spectra to explain the data, 
(b) The corresponding evolution of the radio synchrotron and 
decay 7 -ray flux. The energy ranges of flux integration for the 
radio and 7 -ray photons are 70 MeV—20 GeV and 70 MHz—10 GHz 
respectively. 

pressure in the shell, and the cool shell is supported by 
the non-thermal pressure and stops contracting. In this 
model, the non-thermal pressure is dominated by the 
(transverse) i3-field with Pcr/Pb± ^ 0 . 01 . Compres¬ 
sion of the gas, CRs and R-field then halts in the shell. 
The pressure ratio Pth/Ps (plasma-/3) before and after 
the collapse of the cold shell for this model is < 100 and 
< 1 respectively inside the shell. Since the plasma-/? is 
proportional to nT/B\ and B± reacts to the collapse 
roughly as y/n, its evolution depends mainly on the tem¬ 
perature change only. Since radiative losses of CRs in 
the interesting energy bands are also slow compared to 


the dynamical timescale, the non-thermal emission fluxes 
which mainly originate from the dense shell become sta¬ 
ble after the transition. We note however that we do not 
consider the escape of CRs from the dense shell in our 
model which can result in a decay of the non-thermal 
fluxes with time. However, the emission here are domi¬ 
nated by protons and electrons in the GeV energy range 
embedded in a compressed i?_L of a few 100 /rC whose 
Larmor radii should be very small, and hence we may 
expect them to diffuse away from the shell only slowly 
as well. On the other hand, the continued production 
and accumulation of secondary e“'"e“ from hadronic in¬ 
teractions by the trapped CR protons in the dense shell 
gives rise to a gradual increase of the secondary syn¬ 
chrotron emission after the transition. The total frac¬ 
tion of the initial cloud shock energy attributed to CR 
re-acceleration is < 1 % at the end of simulation in this 
model. 

In the same figure we have overlaid onto our model the 


radio telescopes ( 

Castelletti et al. 

0 

0 

CN 

Sun et aL||2011|) 

and 7 -rav instruments including t'ermi 

jAF (Ackermann 


et al. 20131 and AGILE (Cardillo et al. 20141. Wehndan 


overall satisf actory agreement ot the broadband spectral 
shape. From Castelletti et al. (20071, the radio spectral 


index a = —0.37 ± 0.02, and the range of index from 
our model after the flux has reached the maximum at 
tc > 2 ttr is a « —0.38 to —0.40, which are reasonably 
consistent with each other or at most only slightly softer 
for the model. The predicted radio spectrum shows a 
spectral softening above ^ 10 GHz due to synchrotron 
loss, and is c onsistent with the 70 GHz data from Planck 
observation (Planck Collaboration et al. 2014), but it 
cannot explain the 30 and 44 GHz points which do not 
seem to agree with a simple extrapolation from the lower- 
frequency data. This apparent ‘excess’ can possibly arise 
from the anomalo us microwave emis sion from small spin¬ 
ning dust grains (Scaife et al. 2007) since W44 lies at a 
low Galactic latitude. In this model, most of the bright 
TT^-decay y-rays and radio synchrotron emission originate 
from the cool dense shell which is a thin compressed re¬ 
gion behind the radiative cloud shock. This is qualita¬ 
tively consistent with the filamentary morphology of the 
shell-like remnant revealed in the radio waveband. And 
as a result, the normalization of the non-thermal spectra 
for this model mainly depends on the compression of the 
cold dense shell, rather than the acceleration efficiency. 
To conform with the observed flux level, an overall nor¬ 
malization factor of 0.2 is applied to the model, which 
can be interpreted as a reasonable filling factor of the 7 - 
ray and radio continuum emitting region over the whole 
Att shell. 

It is interesting to note that, in spite of a similar basic 
scenario, a difference between our model and the ana¬ 
lytical results of UlO can be found in the ratio between 
the secondary and primary synchrotron flux in the ra¬ 
dio waveband. They predict a significant contribution 
of the secondary component in the longer wavelength 
regime. Our model predicts a much less prominent sec¬ 
ondary component, and the primary emission alone can 
explain the observed spectrum reasonably well. We be¬ 
lieve this difference stems from the following. UlO as¬ 
sume the presence of a non-evolving, time-independent 
cool dense shell formed behind a radiative shock which 















































9 



Fig. 4. — Brightening of forbidden lines behind the cloud shock 
for the re-acceleration (top panel) and NLDSA model (bottom 
panel). For each line, the emissivities are averaged over the volume 
of the shocked cloud material at each time snapshot. The lines are 
the weighted total of their doublets, if applicable. 

adiabatically compresses all injected primary particles 
once they are re-accelerated at the shock. Re-accelerated 
protons are continuously injected into such a cool dense 
shell and stay there to produce secondary at a fast 
rate (tpp oc up to the present day. Our model, on 

the other hand, treats the evolution of both the primary 
and secondary populations behind the shock fully 

coupled to the hydrodynamics, including the transition 
of the cloud shock into its radiative stage with a time- 
dependent formation of the cool dense shell, which cer¬ 
tainly has an important impact on the production and ac¬ 
cumulation of the secondary leptons through post-shock 
hadronic interactions. 

3.1.2. ‘Forbidden’ Line Emission 

Another interesting order-of-magnitude check we can 
perform is to compare the predicted emissivities of ‘for¬ 
bidden’ optical and IR collisionally excited lines behind 
the radiative shock to observation. Figure |4 shows the 
time evolution of the surface brightnesses ra relatively 
luminous lines included in our calculation of the cooling 
function, up to tc = 2.0 when the computed radio and 
7 -ray flux broadly agree with observations. For the re¬ 
acceleration model, like the y-ray and radio emission, the 
calculated forbidden line emission peak spatially within 
a thin region slightly behind the cloud shock either in¬ 
side the cold dense shell or the adjacent rapidly cooling 


layer, resulting in a filamentary morphology following the 
radiative shock front. 

The re-acceleration model predicts moderately strong 
[O I] AA(6303-b6363) and [O III] AA(4959-b5007) emis¬ 
sion, but observations of W44 in [O I] are not readily 
available in the literature (although very bright 63//m 
[O I] line emission a t the NE rim has been reported by 


Reach & Rho (1996)), and the SNR is highly absorbed, 
s o IU III I emission is unlikely to be detected. 


Mavromatakis et al.| (2003) measured an average ob¬ 
served surface brightness F[s nj « dxlO”^"^ erg cm“^ 
s“^ arcsec”^. With their suggested E(i3 — V) of 3.28, 
the attenuation factor by interstellar absorption is about 


1700 at 6720 A, so the corrected surface brightness be¬ 
comes F[g n] ~ 7x10“^^ erg cm“^ s“^ arcsec”^. Our 
model predicts a [S II] AA(6736+6731) intrinsic surface 
brightness of order < 10“^'^ erg cm“^ s“^ arcsec“^. Con¬ 
sidering the effect of limb brightening from projection, 
this brightness should be further enhanced by a factor of 
a few or an order-of-magnitude, making it fairly consis¬ 
tent with observation. We note that Mavromatakis et al. 
(2003) obtained electron densities below 220 cm“^ from 
the [S II] doublet ratio, while our models produce much 
higher densities because they start with pre-shock densi¬ 
ties of 200 cm“^. This is a problem for any shock model 
in that a high density is needed to account for the flux, 
but a low density is needed to match the [S II] doublet 
ratio. High magnetic and CR pressures help somewhat, 
but some combination of lower pre-shock density, emis¬ 
sion from the photoionization precursor and ionization 
and heating of the cold shell may be needed to match 
both the flux and electron density derived from [S II]. 


3.2. Case of NLDSA with Thermal Injection 


Here in the NLDSA scenario, we consider an efficient 
injection of particles from the thermal pool into the DSA 
process. It is not obvious whether such an efficient injec¬ 
tion can occur at a cloud shock, but we suppose that it 
can happen here so that we can study its observational 
consequences. We fix the so-called injection parameter 
6 nj = 3.7, which means the injection momentnm above 
which the downstream thermal particles are injected is 


Pinj = 6 njPth, 2 , where pth ,2 = y^ 2 mpkT 2 and subscript 
‘ 2 ’ means values taken immediately downstream from 
the subshock. The electron-to-proton number ratio at 
relativistic energies K^p is set at 0.003 to fit the radio- 
to-y-ray flux ratio. Of course, when the shock has slowed 
down to a point where only a very low ion fraction ex¬ 
ists in the photoionization precursor (below a few tens 
of km s“^ as shown in Figure [^b)) and the downstream 
electron temperature is too low for collisional ionization 
to occur at an appreciable rate, densities of protons and 
electrons should become too low for efficient injection to 
occur. In such phases, however, acceleration should be 
very inefficient anyway due to the very small Usk- The 
introduction of a spectral break by ion-neutral collisions 
at low Usk is also effective on quenching efficient non¬ 
linear acceleration. As pointed out by [Maikov et al.| 
(2011), pbr essentially acts as Pmax in the conventional 
NLDSA picture since most of the CR partial pressure is 
contributed by the highest energy particles. The break 
hence causes a large reduction of CR back-reaction and 
shock compression ratio. Therefore, contributions to the 
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non-thermal emission in this phase are relatively unim¬ 
portant. Most of the high-energy particles contributing 
to the non-thermal emission are produced in the earlier 
phase in the evolution of the cloud shock, right after the 
blastwave penetration across the wind-cloud interface. 

We perform the same simulation as the re-acceleration 
case and plot the resulting broadband SED in Figure 
In the results, we find a number of interesting differences 
with the former model: 

The first difference is that both the 7 -ray and radio 
flux do not experience a boost after the transition of 
the shock to radiative phase as drastic as we witness in 
the previous case, only a factor of a few here versus a 
few 100 in the re-acceleration model. In this particular 
model, the instantaneous efficiency of DSA (i.e. frac¬ 
tion of incoming energy flux in the shock rest frame con¬ 
verted into relativistic particles), reaches 70 to 80% max¬ 
imum during the early phase of the cloud shock evolution 
when the shock velocity is still high, and the pressure ra¬ 
tio between the non-thermal components (CR and mag¬ 
netic) and the total (i.e. non-thermal plus thermal), i.e., 
(^’cR + -PBx)/-Ptot, reaches about 0.7 downstream before 
the shock becomes radiative, and close to 1 after the gas 
has cooled down rapidly through radiating away its in¬ 
ternal energy. As expected, the non-thermal emission are 
much stronger than the re-acceleration case in the initial 
phase. However, we find that the high non-thermal pres¬ 
sure can support the radiatively cooling gas and prevent 
it from collapsing into a prominent dense shell. As a re¬ 
sult, the effect from rapid compression is not as strong 
in this scenario. The radio-to- 7 -ray flux ratio here is 
very similar to that in the re-acceleration model since the 
adopted K^p = 0.003 is close to the number ratio between 
the Galactic CR electrons and protons specified in Equa¬ 
tion]^ In this model, the non-thermal pressure is about 
equally shared by the CRs and with Pcr/Pb± ^ 1 
compared to roughly 0.01 in the re-acceleration model 
above. The post-shock plasma-/3 before and after the ra¬ 
diative transition is ^ 1 and > 0.01 respectively. The ab¬ 
solute values are much lower than the re-acceleration case 
due to the much larger CR-amplified R-field and lower 
shocked temperature as a result of efficient DSA in the 
early phase. The total fraction of the initial cloud shock 
energy attributed to CR acceleration is about 33% at the 
end of simulation in this model, much larger than the 
‘test-particle’ re-acceleration model we discussed above. 


The second difference is that the shock becomes radia¬ 
tive at a much earlier time than the previous scenario (at 
tc « 0.5 ttr versus 1.2 ftr)- This can be easily explained 
by the significantly lower post-shock temperature at any 
given time due to a much larger portion of the shock ki¬ 
netic energy being converted into non-thermal particles 
by efficient DSA. As a result, the shock does not have to 
decelerate to a velocity as low as when DSA is inefficient 
to achieve a downstream temperature at which the gas 
can cool down effectively throu gh radiation. A sim ilar re¬ 
sult has also been obtained by Bykov et al. (2013). This 
also means that when the shock became radiative, the 
shock was still strong enough to support NLDSA so that 
at that time Pcr has already dominated the post-shock 
pressure to support the gas from collapsing drastically. 
The different timing of the transition to radiative phase 
is also reflected in the evolution of the forbidden lines; 



loQio E [keV] 


(a) Evolution of Broadband Spectrum 



(b) Evolution of Integrated Flux 

Fig. 5. — The same as FigureJ^but for the NLDSA scenario. An 
overall normalization factor of 0105 is applied to the model spectra 
to explain the data. 


brightening of the lines in the NLDSA model occurs at a 
substantially earlier phase than the re-acceleration model 
(bottom panel in Figure]^. 

The third difference can be found in the radio spec¬ 
trum. The calculated 7 -ray spectrum provides a good fit 
with W44 data, but the spectral index of the radio syn¬ 
chrotron emission is substantially steeper (a « —0.53 to 
—0.54) than the former model and typical observed val¬ 
ues. This model has a much more interesting time evo¬ 
lution of the shock compression ratio as well as the CR- 
driven amplification of R-field in the precursor, which 
are both quickly decreasing with time as the shock de¬ 
celerates in the dense medium; the compression ratio de¬ 
creases from about 8 to below 4 from = 0 to 3 ttri and 
even more dramatically for the amplified downstream 
B-field which drops from almost 1 mC to only slightly 
higher than the compressed ambient field of the cloud. 
This results into the initial gradual decrease of both the 
TT^-decay 7 -rays and synchrotron radio fluxes before the 
transition to the radiative phase. Although the well- 
known NLDSA effect of spectral hardening is witnessed 
for the accelerated electrons in the early stage, quickly 
DSA becomes inefficient and the newly accelerated elec¬ 
trons possess a spectral index steeper than 2 in energy 












11 


space. The formation of a prominent dense cool shell 
which compresses and re-energizes the downstream CRs 
like in the previous scenario fails to happen here as well. 
Consequently the overall electron spectrum is not hard 
enough to reproduce the observed radio index. This is 
even more true if we adopt a larger injection parame¬ 
ter (i.e., less efficient injection) and in such a case even 
the observed 7 -ray spectrum cannot be explained due to 
a weaker B-held amplification and hence lower Pmax for 
the CR protons. 

The secondary-to-primary ratio of radio synchrotron 
emission is also higher throughout the evolution in this 
model. This can be explained by the fact that the pre¬ 
vious model, which does not (re-)accelerate CRs so ef¬ 
ficiently, relies on the formation of a dense cool shell 
later on in the evolution to boost the production rate 
of secondary leptons, while in this scenario protons are 
accelerated much more efficiently from the start to inter¬ 
act with the more strongly compressed post-shock gas. 
For this reason as well as that discussed in the previous 
paragraph, the secondary flux also rises by a much larger 
factor in the re-acceleration scenario through the range 
of time shown. 

An overall normalization factor of 0.05 has to be ap¬ 
plied to this model to reproduce the observed absolute 
fluxes, which is smaller than the former case and may 
not be easily attributable solely to the filling factor of 
the emission region. Since a prominent cool dense shell 
cannot be formed, the brightness profiles (along the di¬ 
rection perpendicular to the shock front) of the emission 
are also more diffuse and less filament-like than the for¬ 
mer case. This is also true for the calculated optical and 
IR forbidden line emission. 


3.3. Beyond Spherical Symmetry and Other Aspects 

So far, we have investigated various interesting physical 
processes and their effects on the non-thermal emission 
at a cloud shock moving in a uniform dense medium, 
and compared them in the context of two particle ac¬ 
celeration scenarios for the emission mechanism. Since 
these discussions are based on models constructed by the 
one-dimensional CR-hydro-NEI code, it is important to 
point out here several expected modifications of our re¬ 
sults when the situation deviates from spherical symme¬ 
try. Although molecular clouds interacting with middle- 
aged SNRs are usually found to be more homogeneous 
than those interacting with young SNRs, discussion of 
these effects is nonetheless of significant interest. 

First of all, the molecular cloud can be clumpy and 
multi-phase. Although molecular cloud formation re¬ 
mains a highly active research field and many uncer¬ 
tainties still exist, one usually expects a lower density 
interclump medium surrounding some higher density gas 
clumps, most probably due to ther mal instability and 


1969 


umpy. 


runaway cooling of the ISM (e.g.. Field et al. 
Koyama fc Inutsuka|2000 ). If the cloud is highly c. 
it is possible that the cloud shock(s) cannot propagate 
in a dense medium for a duration sufficiently long to 
become fully radiative and to produce prominent dense 
cold shells. Also, shock-cloud interactions can happen at 
different times at different clouds such that, at a given 
SNR age, the integrated non-thermal emission becomes a 
superposition from multiple shocked clouds with distinct 
tc, resulting in a change of the overall spectral properties. 


Other complexities like reflected shocks an d MHD tur¬ 
bulen ces from shock-cloud interactions (e.g., Inoue et al.j 


20121 can also happen. 

DsA at shocks in such a highly clumpy medium is ex¬ 
pected to be modified compared to a spherically sym¬ 
metric model, such as a modified CR escape process and 
magnetic field distribution. In the interclump gas with a 
lower density, cloud shocks can propagate at much faster 
speeds than in the dense clumps, so that possibly they 
can accelerate high-energy electrons capable of emitting 
localized non-thermal X-ray em ission, as has been sug¬ 
gested to be the case for W44 (|Uchida et al.||2012|). In 
such a case, depending on the detailed cloud properties, 
7 -rays with a leptonic origin such as inverse-Compton 
scatterings may play a significant role in producing the 
observed high-energy emission. As the next extension 
of the current work, it is fruitful to perform a series of 
one-dimensional calculations to encompass both fast in¬ 
terclump shocks as well as the slower radiative shocks dis¬ 
cussed here, which can be merged to realize a more com¬ 
plete representation of the situation in a clumpy cloud. 
In the further future, we foresee a generalization of our 
code to three-dimensions while maintaining it to be self- 
consistently coupled to DSA and CR transport. 

Finally, 7 -ray emission from “CR-illuminated” clouds 
recen tly detected from the v i cinity of middle-aged S NRs 
(e.g., Uchiyama et al. 2012 Hanabata et al. 2014) has 
not been discussed here. A self-consistent model of such 


emission together with the shell com ponent (see e.g., Lee 
|et al.|[2TO^ [Ellison fc Bykov]|201H for such models of 
young SNRs) is crucial for providing further constraints 
on the particle acceleration mechanism, such as the de¬ 
tailed processes of escape and propagation of CRs from 


the s hock (e.g., Aharonian & Atoyan 1996 Ohira et al. 


2010 1 . Anisotropic diffusion is expected to play an im¬ 


portant role in determining the 7 -ray flux and morphol¬ 


ogy from these illuminated clouds (e.g., Maikov et al. 


2013 Nava & Gabici 2013). The inclusion of this com¬ 


ponent tor middle-aged SINRs is another important next 
step for the work described in this paper. 


4. SUMMARY 

We explored the possibility of explaining the bright 
non-thermal radio and GeV emission from middle-aged 
SNRs through a scenario with re-acceleration of Galac¬ 
tic CRs or NLDSA of thermally injected particles by a 
fast radiative cloud shock, under the assumption that 
a sufficiently strong magnetic turbulence is present in 
the cloud to support a Bohm-like diffusion. Using the 
CR-hydro-NEI code, we followed the hydrodynamics of 
a shock propagating in a gas cloud with a density typical 
of those found in GeV-bright middle-aged SNRs, and at 
the same time calculated the time evolution of the (re- 
) accelerated CRs in the downstream and the generation 
of associated broadband non-thermal emission. 

In the re-acceleration scenario, we found that the tran¬ 
sition of the cloud shock into the radiative phase, and the 
accompanying hydrodynamical effects on the post-shock 
conditions, play crucial roles in producing the bright non- 
thermal emission from the shell of these objects. The 
predicted general properties of the non-thermal emis¬ 
sion agree well with observations and are illustrated by a 
comparison with the broadband spectrum of W44. The 
model produces strong oxygen lines such as [O I] and [O 





































12 


III], and the predicted surface brightness of [S II] agrees 
broadly with the observed value for W44, offering further 
support to the model. 

As for the NLDSA scenario with ‘thermal leakage’ in¬ 
jection, a number of difficulties are found to explain ob¬ 
servations. For example, we find that although the calcu¬ 
lated 7 -ray emission can reproduce the observed spectral 
characteristics quite well, the model fails to explain the 
radio continuum emission, with a spectral index too steep 
to reproduce typical observed values (e.g., W44 as well 
as others like IC 443 with similar indices), even though 
the model considered invokes a highly efficient injection 
for DSA, significant shock modification and hence spec¬ 
tral hardening. Models with less efficient injection are 
unable to explain observations due to even steeper radio 
indices and lower 7 -ray cut-off energies. 

A number of other interesting differences between the 
behaviors of the two models are also discussed. For 
example, in the NLDSA model, the high non-thermal 
pressure in the post-shock gas is found to be an ef¬ 
fective counter force to the rapid formation of a dense 
cool shell during the radiative cooling of the gas, there¬ 
fore rapid compression plays a smaller role in generating 
the observed broadband emission compared to the re¬ 
acceleration case. The transition to the radiative phase 
also happens substantially earlier and is reflected by the 
time evolution of the brightness of forbidden lines. 

We have limited our discussion in this paper to the 
study of general properties of non-thermal emission from 
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